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Introduction 


This is an attempt to construct simple models of the electric field of 
the magnet osphe re , based on a combination of observation and plausible 
theory* 


The origin of the electric field is assumed to be in interplanetary 
space , "where a field - v X B is created by the motion of the solar 
■wind with velocity v relative to the earth* An open magnetosphere is 
assumed, in which some interplanetary magnetic field lines reach the top 
of the atmosphere in the vicinity of the earth's magnetic poles (the 
name polar caps will be used for the regions to which such field lines 
are connected) , conveying with them the interplanetary electric field* 
Observations indicate that this field is quite variable [ Cauffman and 
Gurnett , 1972 ] but that on the average it points from dawn to dusk 
within a circular polar cap centered on the magnetic pole^e.g* Heppner 

1272 ] . 

In this work, the radius ? of the polar cap will he taken as a/4 , 
where a is the radius of the earth. Equipotential lines of the average 
field inside the polar caps - assuming that it is constant and in the 
dawn-dusk direction - are then aligned with the noon-midnight direction 
(Figure l). Outside the polar cap the conducting ionosphere will extend 
the equipotentials to closed magnetic field lines, creating a two-celled 
pattern there with corresponding ionospheric currents. 

* 

Detailed models for this field will be derived in several stages, in 
all of which the conductivity along field lines will be assumed to be high 
enough to ensure the vanishing of E*B everywhere except in the ionosphere. 
At first the rotation of the earth is ignored completely and a simple model 
is constructed which fits the observed properties listed earlier* Next 
the rotation of the earth is taken into account, but the field is assumed 
to be that of a magnetic dipole rotating around its symmetry axis; this 
allows the concept of the electric potential to be retained, which permits 
the derivation of interesting properties, including the use of a conjugate 
potential u which paces the drift of charged particles in the field* 
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Finally, the general case, involving asymmetrical rotation, is briefly- 
discussed. 

Simple Model Without Rotation 


If there exists no time dependence the electric field may be derived 

from a scalar potential 4> and the assumption that in the magnetosphere 
E*B = 0 then reduces to 


<j> '= <j) ( * , p ) (1) 

-where ( , P ) are Euler potentials of the magnetic field £ Stern, 1970 

and references cited there ] . Assuming the magnetic field to be that of 

a magnetic dipole and denoting by a the radius of the earth and by gx 
the dipole term in the expansion of the scalar potential ( g°<v -5 • 10 5 MKS) 

y = a g° (a/r) 2 cosg (2) 


one finds 


U. = a gx (a/r) sin £ 0 

? - 


(3) 


The solution for $ could nov "be derived as follows# Assuming an exter- 
nally generated electric field above the polar cap ionosphere ( < <<) 

which is constant and directed from dawn to dusk, one can solve for 0 
in the ionosphere (where it depends not only on ot and but on a third 
coordinate as well) both for the polar cap ionosphere and for the fringing 
values of 0 which spill out over the edges of the polar caps into the 
sub-polar ionosphere# The values of 0 in the latter case will then deter- 
mine the potential on closed field lines: due to the electric field on such 
lines the magnet ospheric plasma will experience a drift motion, leading to 
general convective motion in the magnetosphere# 
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This leads to a fairly complicated boundary value problem, which will 
not be treated here (the input field used is already very simplified, anyway). 
Instead a form of <fi will be assumed which appears to be in agreement with 
observations. This form is 

sin ( JVa) (4-a) 

sin ( jVa) * Z-oSq ^ ~' b ' 1 

with k an appropriately chosen constant. 


* - ■ K < < V ^> k/2 


To derive some intuitive insight into this, let the polar cap he 

approximated by a plane surface, tangential to the magnetic pole. In this 

plane, let ( ? , ) he plane polar coordinates and let ? = P he 

v v o 

the boundary of the polar cap. Since f = a sine , one finds 

? / ? „ - ( <* / <* o ) V* (5) 

Thus for $ < 2 
J o 

t - - i (P/f 0 ) Bln'f = - 4> 0 v/! 0 (6) 

which gives a constant electric field E q = 4^/ J 4 in the y direction. 

For f > P , on the other hand 

4> = - <t> Q (? Q /- P sin Y ( 7 ) 


giving a field which drops off approximately as 
of k give the familiar two- celled struo 


- (k+ l) 


4 Figure 1 )j 
cturef: for k = 1 


All values 
this resembles 


the well known two dimensional solution satisfying Laplace*s equation 

for o > 9 and continuing a constant field which exists in ? ^ P 
o v ” o 

£ e.go Panofsky and Phillip s , 1955 , section 4-9 ] , while higher values 

of k progressively compress the equipotentials towards the circle P = P 

- v o 
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The choice k = 2 leads to a constant electric field in the equatorial 
plane ( as assumed by Chen £ lg^O ] ) , but observations from 0G0 6 
on the profile of |E( suggest k = 4 . To fit observations, we shall also 

take 


E = 0.02 volt /meter 

o 



In a dipole field is given by ( 3 ) and therefore (4-b) leads in 
the equatorial plane to 

<j> = - <f> ( oL /a g° ) ^ (r/a) sin ¥ 


- 


( d.Ja. g° 


: >v. 


(r/a) 


(k-2 ) 


(8) 


If k = 2 this field is uniform and equipotentials are all parallel to the 
noon-midnight direction. For higher values of k (as are assumed here) the 
equipotentials still tend t^ follow the noon-midnight direction, but they 
curve towards the noon-midnight line, their closest approach to that line 
oc curing as they cross the dawn-dusk axis. A schematic description of these 
lines for k = 2 is given in Figure 2 , while Figure 3 shows them for k = 4 
Axisymmetric Rotati on 


If we assume that the earth and the material threaded by closed field 
lines are undergoing steady rotation with angular velocity around 

the symmetry axis, the electric field can still be represented by a scalar 
potential, since ^/’"Dt — 0 * Assume then that equations (4) represent 
the electric field E* in the rotating frame. Then the electric field 
observed in a non- rotating frame - e.g. by a charged particle entering 
the magnetosphere - is 

E = E* - v X B (9) 

■where v in this case represents the velocity of co- rot at ion® One then 
finds ^ given o( = •>( (r ; 0 ) , ^ - a'f 
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v X B = a u- r 


sin£ i * (Vc<(r,0)* ) 


= V (a u> U) 


( 10 ) 


Hence on closed field lines 

k, 

(f) = - ( o< Q /oi ) / 2 sin( p/a) + atO^ 


( 11 ) 


As far as the rotation of the polar caps is concerned, it will tend 
to twist open field lines into helical shape. The total amount of twisting 
experienced by such lines does not seem to be large, since the length of 
the geomagnetic tail is apparently only of the order of 500 a . An open 
polar field line sweeping past the earth with the velocity of the solar 
wind spends only about 5 hours in contact with the polar cap, from the 
time it merges with a closed field line to the time it again breaks connec- 
tion with the earth's field, and during that time the earth will only rotate 
through 1*5° . 


Particle Motion 


low 

Consider a particle of very low energy, so^that its drift velocity in 
the magnetic field is essentially the electric drift velocity 

v a = B~ 2 ( E X B ) (12) 

Then in the present case 

v = 0 (13) 

showing that such particles stay on surfaces of constant ^ . In particular, 

in the fields of the present model the electric potential in the equatorial 
plane is 

Js k 

<f> = - f Q ( °<o/ a g i ) 2 ( r / & ) 2 sin ^ + a 2 g° (a/r) (14) 

A low energy particle moving in the equatorial plane will conserve the 
above quantity. For very small r , the second term dominates and equi- 
potential lines in the equatorial plane approximate closed circles around 
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around the origin. For large values of r the second term is negligible 
and the first one determines the structure of the equipotentials; in parti- 
cular, for k ^ 2 we obtain a bunch of equipotentials generally aligned 
with the noon-midnight direction, pinched together near the earth if k > 2 . 

These two regimes are separated by the "last closed equipotential" L , 
which also extends to infinity and has the form schematically shown in 
Figure 2 . At the point marked P in the figure, the equipotential surface 
crosses itself, and since the direction of v4> at P is thus undefined, 
this gradient must vanish there. From (lh) , the vanishing of ''if/'i ' f 
implies 

sin ''f = - 1 ( 15 ) 

while the radial component gives 

x k 

(k 4^/2) ( °^ o /a g°) 2 (r/a) 2 sin ''f = - a 2 g° U} (a/r) (l6) 


clearly requires the positive sign in equation ( 15 ), i.e. P is on 
the dusk side. Using numbers introduced earlier (including k = 4) gives 
for P 

r = 6 a 

which is a little high for the distance to the plasmapause bulge - which 
is the usual interpretation of the point P [ Brice . 1967 ] - but has the 
proper order of magnitude. 

The Conjugate Potential 


The Euler potentials ) are not unique and other potentials 

ot 1 ( °^ , ? ) , P 1 ( * , p ) will also describe the field, provided 

^ P')/ 7) ( ) a 1 (17) 

Furthermore, given a well-behaved function oL' ( d , p ) , some ) 

can be generally found for which ( 17 } holds. In particular, if the electric 
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potential <fi («<, ^ ) is regarded as art Euler potential, a conjugate potent i; 
u( **, f* ) will exist such that 

B = V<£ x V u (18) 

It is easy to calculate the rate at which u changes for a particle 
moving with the electric drift velocity v ^ : 


du/dt = v - Vu - - B" 2 ( x B)» Vu « 1 ( 19 ) 

Thus the rate is constant at all places and for all times, and a 
collection of particles which started with a given constant u at t = 0 
will share a surface of constant u at all subsequent times. In other 
words, just as surfaces of constant <f> guide the drift of particles along 
them, surfaces of constant u pace the rate at which such particles 
advance . 

It may appear somewhat surprising that the same function u('X,P ) 
paces the electric drift motion of all particles, regardless of their 
equatorial pitch angle . Consider for example two such particles 

starting at t = 0 from the same field line, one with * 0 = 71 / 2 and 
the other one with a relatively small A . It is not a-priori clear 
here why the two always stay on the same field line: the result is however 
true and the reason is that the electric drift velocity is a magnetic 
field line velocity [ Newcomb , 1956 i Stem , 1966 ; Stern , I 97 O ] . 

If gradient drifts and curvature drifts are included in the calculation 
of the motion, u is no longer conserved and particles with different 
equatorial pitch angles tend to drift onto different field lines. This 
type of motion, for the case of constant equatorial electric field 
(k = 2 ) has been thoroughly investigated by Chen [ 1970 ] 
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In the equatorial plane the two types of surface degenerate into two 
families of lines (generally not orthogonal) - one of them marking the 
trajectories of particles, the other one the successive positions of a 
"front" of advancing particles. 

In order to derive u it is best to consider it as a function not of 
( ot , ^ ) but of ( d. , ^ ) and of . Substituting this dependence in 

(l8) gives, in a few steps 

= - 1 (20) 

* ^ f 

From this u is readily derived by integration, provided 
can be expressed in terms of ( $ , oC ) . 

Because is given by different expressions inside and outside the 

polar cap, the same distinction must be made when u is derived. Suppose 
that inside the polar cap rotation can be neglected and ^ is given by 

<p = -A(o£/o6)^sin v f - au)o<. (21 ) 

This reduces to a constant field in the y direction if plane polar 
coordinates can be used with the approximation = ( * / ^ Q ) 2 * 

Then ( § - a 'f ) 

= - ( <P o /& ) ( Ot / 0(^)2 COS f 

= - £ (<< o< 0 a 2 ) - (4* + ^ Q ) 2 ] 2 

= -( vot - q 2 )2 (22) 

From (20) then 




u 


_ A A 

- q ) 2 a** ® - ( 2 a/<j£> ) ( oi oC )2 cos 'f (23) 


In the plane approximation which is used here, u is a multiple of x , 

which agrees vith the advance of a front of particles given E = E y , 
/s 0 — 

B = - B q z o Rotation can he taken into account by replacing in 

the second term of (2l) vith , hut the result then is no longer 

intuitively simple. 

Outside the polar cap 


k 

<P = - d ( c* )2 si nt - au)o<; (24) 

giving 

r* 

u ~ [ ^ + atO^) 2 j 2 dot (25) 

This must he integrated numerically along a line of constant potential $> 
The only difficulty here occurs near the intersection hetveen the equipoten- 
tial and the y axis, where cos 'f = 0 , causing the denominator - / ^<£/'c>p 
to vanish. However, a transfomation to a new variable ^ = f - ti / z 
shows that the integrand behaves as V in this vicinity and the 

resulting integral therefore converges. 

The function u thus obtained can he conveniently represented in 
the equatorial plane of the dipole, as was done in Figure 3 • In that 
figure curves of constant are drawn solid, while broken lines repre - 

sent lines of constant u « Distances are given in normalized units of the 
order of the dimensions of the plasmapause £ Stern , 1974 , equations 13 
and 24 J . Unit distance here is approximately 0.412 E^ 1 ^ earth 

radii, where E Q = j> Q / ? Q and where « a/4 is assumed, while 

the line u = 0 is at a distance R = 2 , parallel to the dawn-dusk line. 
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One can view the line u = 0 as the starting position of a group of 
very low-energy particles drifting earthward from the geomagnetic tail 
with the E x B drift. Other lines of constni u then mark the posi- 
tion of the field lines to which these particles are attached at selected 
later times „ 

The values of u on these times are proportional to the times required 
to reach these positions, by virtue of equation ( 19 )* It will be noted 
that in the figure the elapsed time between consecutive curves almost 
doubles at each step, because the magnitude e/b of the drift velocity 
increases rapidly as the earth is approached, due to the increase in B . 
The line u = 1 , in the normalized units used here, corresponds to a 

time of. l/2 JT days or about h hours and it will be noted that this 
is the order of the time required here by the particles to reach the 
vicinity of the plasmapauseo 

The method illustrated here can be adapted to the motion of particles 
of arbitrary energy and pitch angle £ Chen and Stem , 1974- ^ ♦ The 

details and qualitative properties of this, however, are beyond the scope 
of this work and will be described in another paper. 
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Asymmetrical Rotation 


If ~DB/<)t does not vanish E can no longer be described by alone. 
Assume as before that the electric field at the top of the polar ionosphere 
is given by - 7^ , with <P the same as in equations (4) . The electric 

field in the magnetosphere is then still given by (9) , with v the velocity 
of the local medium, but the added term - (v x B) is no longer curl-free. 

Suppose that on all closed field lines v corresponds to rigid 
co-rotation, satisfying 

v = uj X r (2 6) 

with u> a constant vector. If 


B = VxA V <* A = 0 (27) 

then it may be shown that 

v x B = V (a«v) + V x (v x A) (28) 

The proof starts from the identities 

vxB = V(v.A) - A*(7xv) - A» 7 v - v» Va (29) 

0 = V x ( v x A ) - v(v«A) •+■ A(v»v) - A° 7 v t v.7A (30) 

When these are added up, the last terms cancel. Also, from (26) 

v = - j V X r 2 c0 (31) 


so that V * v vanishes. Furthermore, by explicit calculation in cartesian 

coordinates with one axis aligned with W 

= 2 


V x v 


(32) 
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A«7v = U X A (33) 

This shows that all terms except those in equation (28) cancel. In axisymmet- 
rical poloidal fields one may choose Euler potentials o(. (r, 0 ) and P = af 
and take 

A - 

Equation (10) is then recovered as a special case. If the earth's field 
is appro ximat ed by a dipole tilted at a small angle relative to the rotation 
axis, this dipole may be resolved into components parallel and perpen- 
dicular to that axis, with corresponding contributions Ai and A 2 to 
the vector potential. Then contributes only to the second term in (28); 

if A 2 is small, it may be neglected, leading again to the electric field 
of a dipole rotating around its symmetry axis. A scalar potential which 
approximately gives the electric field can then be derived and it resembles 
the one given in (ll), except that in the first term there oL refers to 
the entire dipole field while in the second term only the value of oL 
associated with the aligned dipole component is used. 

As a final note, it may be pointed out that a general asymmetric 
magnetosphere may be expected not to undergo purely rigid rotation, since 
its boundary is fixed in the frame of the solar wind. The external field 
of the earth then undergoes a time-dependent change even in the co-rotating 
frame and the actual electric effects are more difficult to analyze. 
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CAPTIONS TO FIGURES 

Figure 1 - The two- celled structure of equipotential lines in the polar 

cap, viewed from above. 

Figure 2 - Schematic view of equipotentials in the dipole equatorial 

plane, for k = 2 . 

Figure 3 - Equipotential lines (solid) and lines of constant u (broken) 

in the equatorial plane of the dipole. The value of k is 
taken to be equal to if . 
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